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ABSTRACT 



We have obtained spectrophotometric observations of 41 anticenter planetary 
nebulae (PNe) located in the disk of the Milky Way. Electron temperatures 
and densities, as well as chemical abundances for He, N, O, Ne, S, CI, and Ar 
were determined. Incorporating these results into our existing database of PN 
abundances yielded a sample of 124 well-observed objects with homogeneously- 
determined abundances extending from 0.9-21 kpc in galactocentric distance. 
We performed a detailed regression analysis which accounted for uncertainties in 
both oxygen abundances and radial distances in order to establish the metallicity 
gradient across the disk to be: 12 + log(0/H) = (9.09 ±. 05) - (0.058 ±.006) x Rg, 
with Rg in kpc. While we see some evidence that the gradient steepens at large 
galactocentric distances, more objects toward the anticenter need to be observed 
in order to confidently establish the true form of the metallicity gradient. We find 
no compelling evidence that the gradient differs between Peimbert Types I and II, 
nor is oxygen abundance related to the vertical distance from the galactic plane. 
Our gradient agrees well with an alogous results for H II regions but is steeper 
than the one recently published by lStanghellini fc Haywood I (120101 ) over a similar 
range in galactocentric distance. A second analysis using PN distances from a 
different source implied a flatter gradient, and we suggest that we have reached 
a confusion limit which can only be resolved with greatly improved distance 
measurements and an understanding of the natural scatter in oxygen abundances. 
Finally, a consideration of recently published chemical evolution models of the 
Galactic disk suggests that reconciling the current range in published oxygen 
gradients is necessary for adequately constraining parameters such as the surface 
density threshold for star formation and the characteristic timescale for disk 
formation. 



^Partially based on observations obtained with the Apache Point Observatory 3.5-nieter telescope, which 
is owned and operated by the Astrophysical Research Consortium. 

^Visiting Astronomer, Cerro Tololo Interamerican Observatory, National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under 
cooperative agreement with the National Science Foundation. 

•^Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under cooperative 
agreement with the National Science Foundation. 

*Now at Department of Astronomy & Astrophysics, University of Michigan, Ann Arbor, MI 48109, USA 
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Subject headings: ISM: abundances, abundances, planetary nebulae: general, 
stars: evolution, galaxies: evolution 



1. Introduction 



The current metallicity of the ISM at any point in the disk of a galaxy serves as an 
endpoint indicator of all the element synthesis processes that have ever occurred in material 
currently at that location up to the present. Therefore, computing a successful one-zone 
chemical evolution model of an arbitrary location using the observed metallicity as a con- 
straint helps us to accurately reconstruct the element enrichment history as well as infer the 
stellar rates of element production associated with the region. 

A chemical evolution model of the entire disk is essentially a linked set of one-zone models 
organized in order of increasing galactocentric distance, while one of the major constraints 
on such a model is the observed behavior of metallicity with radial distance, a behavior 
which is usually found to be continuous, monotonic, arguably linear, and negative., i.e., the 
abundance gradient. Thus models which successfully predict the abundance gradient help us 
understand: (1) the local processes referred to above plus the distribution of matter across 
the disk which fuels them; and (2) the metallicity and rate of infall of halo material which 
formed the disk at each point out from the nucleus. A gen eral review of galac t ic abu ndance 
gradients in disk galaxies and ellipticals can be found in iHenry fc Worthey I (119991 ). while 
many of the papers cited below serve as valuable updates to this article. 

Abundance gradients in the Milky Way disk are commonly tracked by measuring the 
abundance of metals such as O, Ne, S, Ar, and Fe as a function of galactocentric dis- 
tance. The first four elements are m ost readily (but not exc l usively) observed in emis- 



sion line objects s uch as H II regions (iDeharveng et al. 



planetary nebulae (IHenry et al. II2004I : iPerinotto &: Morbidelli 



2000 



Rudolph et al. 112006) and 



2006 



Maciel fc Costa 



2008 



Stanghellini fc Haywood 201o[). Oxygen, the focus of this paper, can also be usefully 



measured in B dwarf stars (ISmartt fc Rolleston Ill997l) . while O and Fe can be measured 
in Cepheid variable stars ( Andrievsky et al. I 2004 ). Most researchers find that the oxy- 
gen abundance gradient is linear in the log-normal plane of O versus galactocentric dis- 
tance and generally range s between -0.01 to -0.07 dex kpc~\ although a few authors, e.g., 
Vilchez &: Esteban I (119961 ) ; iMaciel &: Costa I ( l2008l ) , claim that the gradient may flatten out 
at large radial distances from the Galactic center. Other alpha e lements such as Ne, S, and 
Ar follow similar trends ( iHenry et al. Il2004j : iMaciel et al. II2006I ). 
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Planetary nebulae (PNe) are doubly useful probes of chemical composition, as they 
provide information about low- and intermediate-mass stellar nucleosynthesis affecting He, 
C, N, as well as archiving abundances of heavier elements like O, Ne, Ar, CI, and S at the 
time when the progenitor star formecll]. Here we assume that the progenitor abundance 
is equivalent to the ISM value at the time of star formation. For this second group of 
elements, FN abundances measured at various galactocentric distances have been used fre- 
quently t o determine the abundance gradient in the disks of the Milky Way and other spiral 



Stanghellini fc Haywood 



galaxies flHenrv et al. 1120041: iMaciel et al. Il200a iMagrini et al. 1120071 : iBresohn et al. 112010 



20101 



Recently, FN abundances have been used by Maciel and collaborators [see lMaciel fc Costa 



(120081 ) and references cited therein for a complete summary] to study the evolution of the 
gradient of oxygen and other elements. This particular application follows from the idea 
that since progenitor stars range in mass between 0.08 and 8 M© and therefore collectively 
probe the interstellar medium at different times over billions of years, abundance and radial 
distance can now be coupled with age to infer slope changes with time. As a result, the Ma- 
ciel team claims that the oxyge n gradient and the gradients of N e, S, and Ar have flattened 
with time. In contrast, however, IStanghellini fc Haywood I (120101 . SIO) infer from their study 
that the gradient in the Milky Way disk has steepened with time. Thi s confusion carries 
over into the theoretical realm, where the che mical evolution models by iHou et al. I (|2000[ ) 



predict a temporal flattening, while those by IChiappini et al. I ( 1200 ll ) predict a steepening 
with time. 

In the present paper, we present new spectrophotometric observations of 41 FNe located 
toward the Galactic anticenter with galactocentric distances ranging from 0.9 to 21 kpc. 
Some of these objects have no spectra available in the literature, while others have not been 
observed with modern detectors or over a comparably broad spectral range. These new data 
are used to determine the chemical abundances of He, N, O, Ne, S, CI, and Ar in these objects 
in addition to electron temperatures and densities. Incorporating the inform ation on these 



new objects into our previo usly existing abundance database of Galactic FNe (IHenry et al. 



20041 : iMilingo et al. 1 12010[ ). while adopting distances largely from a single source in the 
literature, provides a completely homogeneous sample of 124 FNe. We employ this database 
to study the Galactic chemical gradient in the disk. 

Our study focuses exclusively on the determination of the oxygen abundance gradient 



^ There i s some evidence that the oriKinal progenitor levels of O and Ne c an be altered as a result of CNO 



processing (jPeauignot et al. II2OOOI : IWang fc Liu Il2008t iMilingo et al. ll201Cll ). However, because this idea has 



yet to be firmly established, in this paper we shall assume that PN oxygen levels are unchanged from the 
levels inherited by the progenitor star from the ISM at birth. 
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of the Milky Way disk through the use of PNe, keeping in mind that the gradient inferred 
from PN studies hkely apphes to an uncertain mix of stars born in an earher epoch than 
the gradie nt measured by H I I regions. While our earlier pap ers have looked at ab undance 



gradients (IHenry et al. Il2004f ) and element-to-element ratios (IMilingo et al. II2OIOI ) for sev- 



eral elements in great detail, for the current study we chose instead to carry out a rigorous 
statistical determination of the oxygen gradient since ions of oxygen, the most abundant 
alpha element, are readily observed, and accurate abundance measurements can be carried 
out. 

Our current study possesses two important characteristics which are relevant to the 
confidence which we will place in the outcome. First, our PN sample is completely homo- 
geneous: we have performed all of the observations ourselves using only three instrumental 
systems (KPNO, CTIO, APO), all data have been carefully reduced in a uniform way, and 
the temperatures, densities, and abundances have been computed in a consistent fashion 
using the same code - and the same atomic data and ionization correction factors - through- 
out. Second, our statistical analysis accounts for uncertainties in both the oxygen abundance 
and the galactocentric distance of each object and uses two independent statistical packages 
which have been checked against each other for consistency and accuracy. 

We present our new observations and the derived physical properties of the 41 new 
PNe, including the abundances, in Section 2. Our statistical study of the oxygen gradient is 
described in Section 3, while Section 4 presents our summary and conclusions. In Appendix 
A we compare our oxygen abundances with values from the literature. 



2. Observations and Reductions 

2.1. APO Observations 

Observations of 37 PNe were carried out on the ARC 3.5m telescope at Apache Point 
Observatory, NM, using the Dual Imaging Spectrograph (DIS). Each side of the DIS has an 
E2V 2048 X 1024 pixel CCD. The DIS allows simultaneous observations to be made in the 
blue and red portions of the spectrum, providing the important advantages of identical slit 
placement and identical sky conditions for the entire optical-near-infrared spectrum from 
3600 to 9600 A. 

We observed with a 360" by 2" slit generally oriented E-W. Counts on the CCD were 
binned by 2 pixels (0."8) along the slit direction before readout. The seeing was generally 
1"-1."5, though for this program of extended objects seeing has little impact except possibly 
for the flux calibration star. On the blue side we used the B400 grating giving 1.83 A/pixel, 
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and a resolution of ~7 A. On the red side we used the R300 grating, yielding 2.31 A/pixel, for 
a resolution of ~9 A. The spectroscopic data are noisy between 5400 and 5700 A where the 
red and blue spectra overlap; in addition, the sensitivity of the CCD decreases at A < 3900 
A and > 8500 A. 

Observations were obtained during five runs in 2007 between January and November. 
We reduced the data in the standard fashion with IRAfEI Bias and overscan levels were 
subtracted and a normalized fiat field was applied. We checked the data for response vari- 
ations in the spatial direction along the slit and found that there was no need to apply an 
illumination correction. For each PN we extracted an appropriate spatial portion of the 2D 
spectrum, collapsing it to one dimension, taking particular care to extract the same size 
spatial swath in the blue and the red, since the spatial scales differ between the blue (0.40 
arcsec/pixel) and the red (0.42 arcsec/pixel) CCDs. The wavelength scale was derived from 
observations of arc lamps, and the flux calibration from observations of two or three standard 
stars each night. Table 1 shows our observing log. 

2.2. KPNO and CTIO Observations 

We also report here new results for four additional PNe included in the gradient analysis 
in §3. These objects, observed in 2003 and 2004 with the Kitt Peak National Observatory 
2.1 m telescope or the Cerro Tololo Interamerican Observatory 1.5 m telescope, are listed 
at the bottom of Table 1. See §2 in Milingo et al. (2010) for details on the instrumental 
configurations and observing procedures. These data were also analyzed in the standard 
fashion using IRAF, as outlined above. 

2.3. Line Measurements 

For each PN, the extracted spectra were combined to make a single blue spectrum and 
a single red spectrum. The only exceptions arose when a strong line was saturated in the 
normal-length exposures: [O III] A5007 in the blue and Ha in the red. In those cases, we 
obtained an additional blue and/or red short-exposure spectrum used to measure just the 
line that was saturated in the longer exposures. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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We measured emission-line fluxes from the final combined, flux-calibrated blue and red 
spectra of each PN. These fluxes formed the input for our abundance determinations. We 
used the ELSA package (Johnson & Levitt 2006) for the subsequent plasma diagnostics 
and abundance calculations. The first step in the analysis is to generate a table of line 
intensities that have been corrected for interstellar reddening and for contamination of the 
hydrogen Balmer lines by coinci dent recombina tion l ines of He"*"^. The redd ening curves 
that we employed w ere taken from lSeaton I (1 19791 . UV). ISavage fc Mathis I (119791 visual), and 
Fluks et al. 1 f ll994j . IR). ELSA iterates calculations of the density (from [S II] A6717/A6731) 
and temperature (from [O III] A5007/A4363) until a convergent solution is reached, then 
calculates the amount of contamination to subtract. It also calculates the correct Ha/H/3 
ratio for the converged values of temperature and density. Table 2 lists the observed fluxes 
and corrected intensities, including their uncertainties, along with the calculated value of c 
(the logarithmic reddening parameter), the appropriate value of Ha/H/3, and the observed 
H/3 flux through the spectrograph slit. 



2.4. Abundances 

Along with the [O III] temperature and [S II] density calculated as described above, 
ELSA calculates the [N II] temperature from A6584/A5755. If the required lines are ob- 
served, the [S III] temperature from A9532/A6312 or A9069/A6312, the [S II] temperature 
from AA6717-F6731/AA4068+4076, the [O II] temperature from A3727/A7323 and the [CI III] 
density from A5517/A5537 are also calculated. For most of the ionic abundance determina- 
tions to be described below, either the [O III] or [N II] temperature was used along with the 
[S II] density. 

Ionic abundances for our sample are given in Table 3. The table shows the emission line 
used to calculate each listed value of ionic abundance along with the diagnostic temperature 
used in the calculation. The designation "wm" refers to the weighted mean of the ionic 
abundances, weighted by the observed flux of the line used to calculate it. Only lines flagged 
with an asterisk are included in the weighted mean. At the end of each element's list of ionic 
abundances in Table 3 is the ionization correction factor (ICF). This factor represents an 
empirical attempt to correct the observed ionic abundances for contributions by unobserved 
ionization states. The ICFs themselves take advantage of similarities in ionization potentials 
among different elements, though for some, like sulfur, there are no close similarities, and 
we rely on model predictions. See, e.g., Kwitter & Henry (2001) for a complete discussion 
of the methods used to compute ionic abundances and the ICFs used in our work. 

Electron densities and temperatures for each object are presented in Table 4, while 
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final elemental abundances and selected element-to-element ratios are given in Table 5. The 
statistical uncertainties provided in all cases are the result of careful propagation by ELSA 
of line strength measurement errors, while systematic uncertainties are not included. 

In Appendix A we compare our oxygen abundances with others from the literature. 
We assert that although the abundance agreement is generally very good, there are some 
outliers. 



3. The Galactic Oxygen Gradient 



3.1. The Determination of the Oxygen Gradient 

Table IH] is a compilation of information for our sample of 124 PNe that is used in the 
analysis to follow. For each object identified in column 1, column 2 provides the heliocentric 
distances taken from Cahn et al 19 92 for each of the objects. All but eight of these values 
were available in lCahn et al. I (119921 ). References for the remaining eight objects are provided 
in a table footnote. The corresponding galactocentric distances in column 3 were computed 
using the formula given in a footnote to Table El where we took the sun's galactocentric 



distance t o be 8.5 kpc. Emp 



the LMC, IStanghellini et al. 



oying a distance method that was calibrated using PNe in 



(120081 SSV) provide updated heliocentric distances for 101 of 



our 124 objects. These values and their corresponding galactocentric distances are listed in 
columns 4 and 5. Finally, column 6 lists the final oxygen abundances from Ta bles 5, using 



the cu stomary form 12-|-log(0/H). The following analysis is based largely on the lCahn et al. 



( I1992I ) distances, but we use the SSV results to demonstrate the effects of using different 
distance scales. 

Fig. [T] is a plot of 12+log(0/H) versus galactocentric distance (column 3, Table E]) in 
kpc for the 124 objects in ou r PN sample with w el l-measured dist a nces. The PN sample is a 



compilation of objects from iHenry et al. I (|2004( ). iMilingo et al. I (120 lOf ). plus those objects 
reported upon for the first time in this paper (HKIO). Our sample contains only those 
objects classified as either Type I or Type II, i.e., objects whose location and kinematics 
indicate that t heir progenitors wer e members of the disk population. We use the heliocentric 
distances from ICahn et al. I (119921 ) from which we calculate galactocentric distances, except 
for the following objects (h eliocentric distance sou rces in p arentheses): H 2-18, He2-48, He2- 



6853, NGC 7293 (iHarris 



St3-1 (ITajitsu &: Tamura 



55 (iMaciel 1119841) : IC 418 (jGuzman et al. 1120091): M3-15 (IZhang 1 119951): N GC 6720 NGC 



20061) : H3-75 (lAmnuel et al. 1119841 ): K3-64, K3-93 ( iJaskot II2OO8I ): 



19981 ). 



Our primary goals are: (1) to subject the data to careful statistical tests in order to 
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determine the nature of the oxygen abundance gradient in the Galactic disk as measured 
by PNe; and (2) to compare our result with analogous studies using other types of probes. 
In order to be included in the gradient analysis, we required that a PN have good determi- 
nations of distance, plus good Te[0 III], Te[N II], and Ne[S II], which means reliable flux 
measurements of A4363, A5755, AA6717/6731. Combining the 41 PNe whose measurements 
are presented here with additional disk PNe already in our database, 124 objects qualified 
for inclusion in the sample on which we performed a thorough linear regression analysis. 

Our first step was to gauge the strength of a correlation between O and Rg by comput- 
ing a Pea rson correlat i on coe fficient, r, for the entire sample of objects, using the program 
pearsn in iPress et al. I (j2003[ ). The Pearson correlation coefficient is the ratio of the covari- 



ance of two variables-in this case O and Rg - to the product of their individual standard 
deviations. The square of this coefficient, r^, called the coefficient of determination, is a 
measure of the strength of the linear relationship between the two variables. Both r and r^ 
are i nvariant under linear transforrn atio n, and thus un certainties in O and Rg are irrelevant. 
[See iRodgers &: Nicewander I (119881 ) and iHays I (Il98ll ) for a more detailed description of the 
correlation coefficient.] 

For our complete PN sample we found that r=-0.54 and thus r^=0.29. The proper 
interpretation of this result is that 29% of the variability in 12+log(0/H) (henceforth O) is 
accounted for by Rg under the assumption that a linear model exists which describes the 
0-Rg relation. The same statistical program also predicts the probability that this result 
could arise from a completely uncorrelated parent population. Thus, if the null hypothesis is 
that O and Rg are uncorrelated in the parent population, then the probability that a sample 
of 124 objects drawn from such a parent population will have a value of |r| >0.54 is only 
3.3 X 10~^^. In other words, there is a vanishingly small probability that this null hypothesis 
is being falsely rejected (Type I error; Hays 1981), and so the correlation appears to be truly 
nonzero. 

Our next step was to derive a good linear model, i.e ., least squares fit, f or the data in the 
0-Rg plane. We did this by using the program fitexy in lPress et al. I (|2003[ ). which accounts 
for errors in both coordinates: O and Rg. Our initial attempt included the one-sigma 
uncertainties in O taken directly from our published estimates, wh ile for Rg we assumed a 
standard one-sigma uncertainty of ±20% ( IStanghellini et al. II2006I ). These factors resulted 
in a slope, b, of -0.066±.0055 dex/kpc and an intercept, a, of 9.15±.04. Assuming that the 
errors are normally (Gaussian) distributed, this fit has an associated value of 178.0 with 
122 degrees of freedom u (124 sample objects minus 2 determined parameters, a and b), or 
a reduced ixl = X^/^) of 1.46. 



Ideally, we would like the value of xl fo be unity; values greater than this suggest that 
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the random errors have been underestimated (or there is an unidentified factor that results 
in the scatter being greater than what is described by the uncertainties). Meanwhile, values 
less than unity indicate that one has overestimated the errors. In the spirit of Goldilocks, 
when ^ 1, then the assumed uncertainties satisfactorily explain the scatter of the data 
points around the model. A method of quantifying this idea is to compute the probability 
that a value of xl greater than or equal to the one determined from the data could result 
from randomly sampling the p arent population described by the model and for the same 
number of degrees of freedom ( iBevington &: Robinson II2003I ). This is done by integrating 
the probability function from xl to infinity to obtain q^2, often referred to as a goodness- 
of-fit parameter. A relatively large q^2 suggests that within the confines of the one-sigma 
errors, xl determined from the data has a reasonable chance of occurring and the model is 
likely adequate in characterizing the data. Likewise, a relatively small value means that the 
errors cannot explain xl and that the model is likel y inadequate . The threshold value for 
q^2 in statistical studies is usually taken to be 0.05 (IHays Ill98ll ). In our current case, q^2 
was determined by the program fitexy to be 0.00074, and so we deemed this particular trial 
model unsatisfactory. 

Before continuing further with our analysis, we decided to write our own program 
for computing line ar least-square fits in order to verify the accuracy of the routines in 



Press et al. I (120031 ) as well as to understand the methods better. Therefore, we wrote a 



program in Mathematica to generate a least-squares fit of data with errors in both c oordi 
nates. Spec ifically, we so lved iteratively the exact "l east-squares cubic" derived by 



York 



fll966l . llQGQf ) as unified bv lYork et al. I f l2004h (see also MacDonald &: Thompson Ill992h . To 
verify our procedures and implementation, and assess the accuracy of the resulting fitting 
paramete r s we successfully checked two test cases: (1) the cla ssic da.ta cas e developed by 



Pearson I (jlQOll ) with er rors assigned to both coordinates bv (lYork 



pre viously considered b y iLybanon I (119841 . Il985l ) ; iPowell &: Macdonald I ( 



966, Table II) and 



1972|);lReedJ (11992) 



and I York et al. I (120041 ): and (2) data for electron temperatu res of ionized hydrogen in the 



Magellanic clouds derived by IVermeii fc van der Hu 



S"*"^ and previously considered by I York et al. 



mm 



St 



(120021 ) from abundances of O"*"^ and 



In attempting to produce a fit to our own data sample of 124 objects, we assumed the 
same one-sigma errors as in the above trial with fitexy. Only six iterations were required 
to determine the slope and intercept to 10 decimal places; the rounded values are b = 
— 0.066±.006 and a = 9.15±.04, a x^ value of 178.0 and a xi of 1.46. The calculated standard 
deviations in the slope and intercept (the "standard errors" ) were computed following the 
guidelines in I Cowan I (119981 Cha p. 9). Since our r esults are in excellent agreement with 



those produced by the routines in lPress et al. 
using the latter formulations. 



(12003! ) ■ we confidently continued our analysis 
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We sought to improve the fit by scahng the uncertainties of both O and Rg upward 
independently, since increasing the standard deviation of either or both coordinates reduces 
the size of correspondingly, the value of q^2 goes up. For each PN the new uncertainty 
(To = ctq X f(0), where ctq is the originally-determined uncertainty in O and f(0) is the 
scaling factor. For Rg, we set ciRg = Rg x f(Rg). Thus, we experimented with different values 
of the scaling factors f(0) and f(Rg) in systematic fashion in order to find reasonable ranges 
in these numbers for satisfactory fits. Values of f(0) ranged between 1 and 2, while values 
of f(Rg) ranged between 0.0 and 0.3. 

The test results for each combination of scaling factors are given in Table [71 where the 
first two columns list the scaling factors, the third column lists the resulting slope of the fit, 
and the fourth column provides the value for the q^2, with the rows ordered top to bottom by 
increasing q^2 . We see that acceptable models correspond to gradient slopes ranging between 
-0.041 and -0.074 dex/kpc. Additional tests not reported here indicated that larger values 
of f(Rg) will push the lower limit downward but it is not likely to go below -0.08. However, 
realistic values of f ( Rg) li kely do not greatly exceed 0.3, based upon the brief discussion in 



Stanghellini et al. I (120061 ). 



To refine ou r model further, we assumed that f(Rg) = 0.2, consistent with the value 



recommended by IStanghellini et al. I ( 120061 ). and then adjusted the value of f(0) until the 
fit explained all of the scatter in O, i.e., x^=l-0. This criterion was met when f(O)=1.40, 
at which point the resulting model had the following parameters: a=9.09±.05 dex, b=- 
0.058d=.006 dex/pc^, x^=121.8, x^=1.00, and q^2 = 0.49. Since the last statistic is well above 
the threshold of 0.05, the model cannot be rejected, and we judge the fit to be acceptable. 
These results for the total PN sample are shown in the first row of Table [HI where the first 
column describes the sample or subsample under consideration, the second column provides 
the number of objects in the sample, and the next seven columns give the y intercept, 
slope, x^i reduced x^j goodness-of-fit parameter, correlation coefficient, and the correlation 
probability factor, respectively. (Subsequent rows provide the statistics for PN subsamples 
or other object types.) 

Thus, we adopt the following analytical behavior of the oxygen abundance with galac- 
tocentric distance (Table [HI columns 3 and 6) as that which best describes the data for the 
total sample presented in Fig. [1] 

12 + log(0/H) = (9.09 ± .05) - (0.058 ± .006) x Rg. (1) 

This regression model is indicated with a solid bold line. At the Sun's distance from the 
Galactic center, 8.5 kpc, our model predicts a value of 8.60±.07 for 12 -|-log(0/H), close to a 



recently-determined solar value of 8.69±.05 by lAsplund et al. I (120091 ). 
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We checked to make sure that our inferred gradient in eq. [T]was not being influenced too 
heavily by the few extreme points at large Rg. To do this, we started with the object with the 
greatest galactocentric distance and excluded PNe one by one, calculating a new value for 
the gradient each time. All resulting gradients obtained in this manner were slightly flatter 
- the maximum slope was -0.051, occurring after eliminating all of the outer six objects. 
Thus, we feel confident that the model in eq. [1] is not being influenced unduly by the few 
objects in our sample located at large radial distances. 

The need to increase the value of ctq by 40% beyond our previously established level 
of uncertainty in order to reach xl = ^ suggests that objects in our PN sample possess 
some natural scatter which is related to real abundance differences among them. Such 
differences could be caused by one or more of the following: (1) age differences of the 
progenitor stars, where the older ones formed out of less metal-rich material than younger 
ones; (2) inhomogeneous mixing of the interstellar material out of which the progenitors 
formed; (3) stellar diffusion, in which stars migrate along the disk from a galactocentric 
radius where they formed to their present location; and/or (4) systematic errors in the 
abundance determination process. It is clear from Table [7] that the value of the gradient is 
sensitive to uncertainties in Rg as well as to the oxygen abundance. However, our adopted 
standard uncertainty of 0.20 for Rg is considered to be quite reasonable (L. Stanghellini, 
private correspondence), and thus we still maintain that at least some natural scatter in the 
oxygen abundances exists. In any case, it is beyond the scope of this paper to further discuss 
the origin of the natural scatter. 

For comparison purposes and to see the effects of using a different set of distances, we 
used the 101 objects in Table [6] for which SSV provided heliocentric measurements in order 
to repeat the regression statistics computation. We first compared directly the two distance 
sets and found that while most objects fell very close to a line of 1:1 correspondence, roughly 
10% of the objects did not. Interestingly, these differences turned out to have a noticeable 
impact on the regression results when performed with the SSV distances. Using the same 
uncertainties that produced the result in equation [1], i.e., f(0)=1.4 and f(Rg)=0.2, the SSV 
distances now imply a slope of -0.042±.004 with x^=1.40. Additional tests involving larger 
(probably unrealistic) values of f(Rg) to drive xt down resulted in a further flattening of the 
slope. We shall comment on this result at the end of this subsection. 

Next, we subdivided our sample by galactocentric distance, Peimbert type, and height 
above the Galactic plane. In the first two cases we computed a least squares fit and correlation 
coefficient for each subgroup to look for distribution differences among them. 

We divided objects into two groups based upon whether their galactocentric distance 
placed them outside or inside the 10 kpc circle. This distance was chosen because it coincides 
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with the position of the sharp discontinuity in [Fe/H] reported by iTwarog et ah I (119971 ) in 
which the mean value of [Fe/H] beyond this point was found to drop by 0.3 dex. We 
then computed the regression statistics separately for the two groups. PNe out to 10 kpc 
have a gradient and correlation coefficient of -0.054±.013 and -0.21, while those beyond 
this distance have values of -0.12±.14 and -0.39, respectively, and their curves are shown in 
Fig. [U In the same figure we also show a quadratic fit to the data which also indicates a 
steepening of the gradient, where the function corresponding to this curve is 12-|-log(0/H) = 
8.81 — 0.014 X Rg — 0.0011 x R^. The quadratic fit predicts a relatively fiat gradient at small 
galactocentric distances but a rather steep one past 10 kpc. Although these results suggest 
that the gradient steepens in the outer disk, we caution that the scatter is broad beyond 
10 kpc, and only additional data at these distances will allow us to reliably determine the true 
behavior of the slope. Inter estingly, some authors present empirical evidence for a flattened 



gradient in the outer disk ( Vflchez fc Esteban 19961; 



20081 : IPedicelh et al. 1120091 ). For example. ICosta et al. 



Costa et al. 2004; Maciel fc Costa 



( 12004| ) use PNe to estimate that the 



oxygen gradient is -0.0 9 dex/kpc between 4 a nd 5 kpc, but the slope becomes fiat at 11 kpc 
and beyond. Likewise, iPedicelli et al. I ( 120091 ). using a sample of 265 Cepheids between 5-17 
kpc from the Galactic center, find the iron gradient to be three times steeper inside the 8 kpc 



circle than outside of it. This disagr eement is a. 



of the chemical evolution m odels bvlFu et al. I (12009) and the models employing a constant 



so ap parent in published model results: All 



star formation efficiency by iMarcon-Uchida et al. I (12010 ) pred ict a steepening gradient in 
the outer disk, while the models by lMarcon-Uchida et al. I (120 lOl ) which assume that the star 
formation efficiency falls off with radial distance, predict a flattening of the gradient. 

Figure [2] shows the distributions by Peimbert type along with the regression lines, while 
the complete statistical information is provided in Table [H The gradient slopes (and uncer- 
tainties) were found to be -0.061±.008 and -0.053±.010 for Type I and Type II, respectively, 
while the correlation coefficients were -0.59 and -0.44, respectively. Thus, within the un- 
certainties the two types appear to be indistinguishable, a conclusion supported visually in 
Fig. [2j In addition, the intercept values are also comparable, and so we are unable to see any 
statistically significant difference in oxygen abundance patterns between the two PN types. 
This conflicts with the claim of a gradient difference between PN types by SIO, who flnd a 
slope of -0.035±.024 dex/kpc for Type I PNe and a slope of -0.023±.005 dex/kpc for Type II 
PNe. These authors continue by relating PN type with progenitor age and conclude that 
the gradient has steepe ned with time. Interestingly, their conclusion conflicts with that of 
Maciel fc Costa I (120081 ). who claim that the gradient has flattened over time. Resolving the 
problem of temporal slope change is extremely important from the standpoint of chemical 
evolution of the MWG, but more than likely the situation is currently being complicated 
by small-number statistics along with real scatter in the oxygen abundances of objects at 
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similar galactocentric distances. 



Figure [3] separates the sample by vertical distance fr om the Gal actic plane, where we 
adopt Z=300 pc, the general disk population scale height (ICox Il2000l ). The gradient slopes 
were found to be -0.038±.007 and -0.033±.009 for objects above and below this height, re- 
spectively, while the correlation coefficients are -0.63 and -0.39, respectively. These numbers, 
along with a visual inspection of the lines in the figure, indicate that the distributions of the 
two subsamples are indistinguishable. 

W e repeated t he above exercise but this time adopted Z=100 pc, the Population I scale 
height fICox Il2000l ). The regression analysis inferred slopes of -0.038±.006 and -0.013±.028 
for PNe above and below 100 pc respectively, while the correlation coefficients are -0.59 and 
-0.09, respectively. Here the differences, particularly in the correlation coefficients, are likely 
the result of the relatively small galactocentric distance range for those objects within 100 
pc of the plane. 

Finally, eq. [1] was used to normalize all oxygen abundances to the same galactocentric 
distance and then these adjusted oxygen abundances were plotted against the corresponding 
z distance perpendicular to the Galactic plane. While we do not show the plot here, we 
can report that no correlation was found. We conclude that within our database there is no 
evidence for a correlation between oxygen abundance and vertical distance from the plane. 

Table [H] also provides comparisons of our derived PN oxygen abundance gradient for 
the disk of the MWG with analogous results for what are likely to be younger populations 
of H II regions and B V stars. The bottom two rows of Table [8] show the results of ap- 
plying our statistical methods to an H II region and a B V star sample which we compiled 
using abundances taken directly from the literatur e . The H II region sample was com piled 



from measuremen t s rep orted in lAfflerbach et al. I (119971 ). IVflchez fc Esteban I (119961 ). and 



Deharveng et al. I (l2000l ). and covers the galactocentric distance range of 3 to 17 kpc. We 
used the abundance uncertainties as quoted in the papers and adopted a standard error in 
galactocentric distance of 2 kpc for an accep tabl e fit. Likewise, the B V star sample was 
compiled from studies by lSmartt et al. I (12001! ) and lRolleston et al. I (12000! ) and extends from 
2-18 kpc. In this case it was necessary to scale the abundance uncertainties by 1.9 in order 
to achieve a suitable fit, while the galactocentric distance errors were taken directly from the 
papers. We present a comparison of the total PN, H II region, and B V star samples listed 
in Table IH along with their corresponding least squares fits, in Figure !4l 

In both the H II region and stellar studies we see strong inverse correlations between 
oxygen abundance and galactocentric distance and good agreement in slope values with 
that of our total PN sample. Note that the correlation coefficients in both cases indicate 
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the presence of a stronger correlation than for our total PN sample, although the gradients 
are consistent. Interestingly, this is despite the fact that there is likely to be a mean age 
difference between our PNe and the H II region and B V samples, as the latter two generally 
represent a younger population than the first group. 

Table [9] presents a broad comparison of our adopted gradient (eq. [1]) for our total sample 
with many other gradients from PN studies along with H II region, Cepheid, and B V star 
studies. In contrast to the gradients listed in Table [HI which we computed directly from the 
published data, the gradients in Table M are those provided by the authors. For the source 
identified in column 1, the second column gives the number of objects in the sample, while 
the next three columns provide the gradient in dex/kpc, the type of objects in the sample, 
and the galactocentric distance range covered, in kpc. 

First, we see a wide range of gradient values among the PN samples listed. In particular, 
the difference between this paper and SIO is interesting, since both samples contain a large 
number of objects and cover a broad range in galactocentric distance. Clearly, the gradients 
inferred by these two studies are statistically different, with the SIO paper proposing a flatter 
gradient than the one we find here. We shall expand on the meaning of this difference in the 
next subsection. In the meantime, we verified that the difference in the linear models derived 
by SIO and us were not attributable to the use of different ox ygen ICFs. We recompu t ed ou r 
abundances using their ICF formula for oxygen, taken from iKingsburgh fc Barlow I ( 1l994l ). 
and found essentially no difference in either slope or intercept. 



The other three PN-derived gradient s bv lHenrv et al. I (12004 H04). ICosta et al. I ( 12004 
C04), and lPerinotto &: Morbidelli I (120061 P06) cover slightly more restricted distance ranges. 
H04 used their own spectrophotometric observations to derive abundances of Types I and 
II PNe in the MWG disk; C04 likewise used their own spectral measurements to derive 
abundanc es of 26 Galactic anticent er disk objects and combined these with PNe from the 
sample of iMaciel Sz Quireza I ( 1l999l ) to infer their gradient; and P06 compiled observations 
from the literature of Type II PNe but reprocessed all line strengths themselves to produce 
a homogeneous abundance set for determining their gradient. The gradients of SIO, H04, 
and P06 appear to be very similar, given the uncertainties associated with each. 

The remai ning gradients in Table [9] include those deri ved from optical and IR studies of 
H II regions bvlRudobh et al. I J2006I. R 06). Cepheids bv IPedicelli et al. I J2009I . P09), and 
B V stars bv lSmartt fc Rolleston I Jl997l . SR97). (Note that P09 report an Fe/H gradient.) 
Gradients in these four cases are consistent with those for PNe, although the uncertainties 
are somewhat larger in the former. 



In summary, we see strong evidence once again for a negative oxygen gradient in the 
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disk of the MWG. However, there is no evidence in our own work that the slope differs 
between PN Types I and II or between objects of different vertical distance from the Galactic 
plane. Because of the broad abundance spread beyond 10 kpc, we cannot confirm with any 
confidence the claim by many others that the gradient changes (either steepens or flattens) 
in the outer regions of the disk. We therefore adopt the simplest hypothesis, to wit, for now 
that the gradient is constant along the disk's entire length. 

The most troubling aspect of our gradient derivation exercise is that both the uncer- 
tainties in the distance scale and the oxygen abundances prevent us from pinning down the 
value of the slope at this time to an accuracy better than about 0.02 dex/kpc. This became 
apparent above when we substituted the SSV distances for the ones used in the original anal- 
ysis, i.e., column 2 of Table El and were forced to add more natural scatter to the oxygen 
abundances in order to reduce the value of xl to an acceptable level. From our complete 
exercise we consider it very likely that the true slope is within the range of -0.04 to -0.06 
dex/kpc, but we cannot refine the number beyond that point. Essentially we have reached 
a confusion limit regardingthe abundance gradient as derived using PNe. Perhaps with the 
launch of the GAIA probqj in 2012 trigonometric distances to many Galactic PNe can be 
measured with much better accuracy, and the uncertainty introduced by the distance scale 
for PNe can be significantly reduced. 



3.2. Gradient Uncertainty and Theoretical Models 

Given the range in values for the slope of the Galactic disk oxygen gradient displayed 
in Table the obvious uncertainties related to an unsettled distance scale, and the exten- 
sive work that is currently going on to reconcile these differences, it is appropriate to ask 
whether continuing our efforts to better define the characteristics of the oxygen abundance 
distribution is likely to offer much in the way of additional improvement in our understand- 
ing of the chemical evolution of the MWG disk. To address this question, we consider the 
model-predicted sensitivity of the Galactic oxygen gradient to a few parameters related to 
galactic chemical evolution. 



Recent detailed chemical evolution models of the Milky Way disk by lMarcon-Uchida et al 



( I2OIOI ) illustrate how the value of the present day gradient as reflected by studies of H II 



regions is influenced by the threshold density for star formation, i.e., the surface density 



■^Thc main goal of the Gaia mission, as stated at the GAIA website, "is to make the largest, most precise 
three-dimensional map of our Galaxy by surveying an unprecedented one per cent of its population of 100 
billion stars." 
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above which star formation occurs but below which it doesn't, and the variation of the star 
formation efficiency, i.e., the ratio of the mass of stars formed to the mass of gas available for 
forming them. Their models predict that between 4-14 kpc in galactocentric distance, when 
the star formation efficiency is held constant and the star formation threshold is reduced 
from 7 to 4 MqPc"^, the predicted oxygen gradient flattens, going from -0.059 to -0.025 
dex/kpc. Note that these gradient values are almost exactly the same as those reported 
in this paper (-0.058) and the one by SIO (-0.023), respectively, for our total samples. At 
the same time, holding the threshold constant at 4 MqPc~^ but allowing the star formation 
efficiency to be a function of galactocentric distance produces little effect. Similar patterns 
are seen over more restricted ranges in Rg. Therefore, these models indicate that improving 
our knowledge of the oxygen gradient in the Galactic disk will allow us to pin down the value 
of at least for one important parameter, the star formation threshold. 

In a similar fashion, detailed chemical evolution models of the Milky Way disk by 



Fu et al. I (120091 ) tested the effects of the infall timescale, the star formation law, and disk 
formation timescale (DFT) on the size of the abundance gradient. In their models the disk 
is assumed to begin forming out of matter infalling from the halo at a time specified by the 
DFT, while the rate of infall decreases by e~^ over a time defined by the infall timescale. At 
the same time, the star formation rate is proportional to some power of the local gas surface 
density and may also be a function of radial distance along the disk. 



Fu et al. I (120091 ) tested the impact of these three parameters on the steepness of the 
present day abundance gradient and found that the DFT had the largest effect, while the 
other two parameters influenced the gradient in a minor way only. For example, when DFT 
is zero (disk formation begins at the moment of Galaxy formation) the gradient ranged in 
value from -0.009 to -0.027 dex/kpc. However, when the DFT value increases radially the 
resulting gradient range steepens to -0.056 to -0.091. The range in values in each case is the 
result of changing the star f ormation law or the infall tim escale. The point here, of course. 



is that as in the case of the iMarcon-Uchida et al. I (120101 ) models we see a predicted range 



in gradient steepness that corresponds roughly to the observed range in the oxygen gradient 
presented in Table M 

Thus, chemical evolution models strongly indicate that some important parameters of 
chemical evolution of the disk, the density threshold for star formation and the timescale for 
disk formation, are expected to have a measurable influence on the value of the abundance 
gradient. Therefore, continuing to reflne our quantitative description of the gradient is 
clearly worthwhile. Further studies should continue to focus on extending the baseline out 
to a radial distance beyond 20 kpc from the Galactic center in order to flrmly establish the 
character of the gradient in terms of its actual value in addition to its shape, since there 



- 18 - 



is observational support for a gradient of constant slope plus evidence for gradients which 
steepen or flatten with increasing galactocentric distance. 

4. Summary/Conclusions 

We have obtained new spectrophotometric measurements of 41 planetary nebulae lo- 
cated in the direction of the MWG anticenter using the facilities of the Apache Point Obser- 
vatory. Density, temperature, and abundance information for each object were subsequently 
derived using the software package ELSA. These results were in turn added to our large 
database of homogeneously determined abundances in order to study the oxygen abundance 
gradient in the Galactic disk. The complete PN sample comprised both Type I and II ob- 
jects, 124 nebulae in all, extending in galactocentric distance from 0.9 to 21 kpc, where all 
observations, data reductions, line measurements, and abundance computations were carried 
out using exactly the same methods and executed by the same personnel throughout. 

The primary focus of this paper (in addition to presenting new spectroscopic data and 
abundances for 41 anticenter PNe) was to carefully study the nature of the oxygen abundance 
distribution across the Galactic disk. To accomplish this, we computed linear regression 
models (least squares fits) of the total sample along with subsets by including uncertainties 
both in the abundances themselves and in the galactocentric distances of the objects. We 
have concluded the following: 

1. Our total sample of PNe is well fit by a linear regression model whose form is 12 + 
log(0/H) = (9.09 ± .05) - (0.058 ± .006) x Rg with a correlation coefficient of -0.54, 
where Rg is galactocentric distance in kpc. In computing this model we accounted for 
uncertainties in both the dependent and independent variables. Related probability 
calculations for testing the strength of these parameters support the validity of our 
results. A parallel analysis using the recent distances for a smaller sample from SSV 
resulted in a oxygen abundance slope of -0.042±.004, i.e., somewhat flatter than the 
above results. 

2. We find clear evidence for natural scatter in oxygen abundance for our sample PN 
population, scatter well beyond what can be explained by observational (statistical) 
uncertainties. 

3. Because of the problems related to distance scale uncertainties and real scatter in 
oxygen abundances, we propose that gradient determination studies involving PNe have 
reached a confusion limit, and that further resolution will only come with improved 
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PN distances and understanding of the natural scatter in oxygen abundances which 
appears to exist. 

4. We see no obvious indications that Type I and Type II PN oxygen abundances are 
characterized by different gradient slopes. Likewise, there seems to be no correlation 
between an object's oxygen abundance and its vertical distance from the Galactic 
plane. 

5. We find some evidence that the oxygen gradient beyond 10 kpc is steeper than it 
is inside this distance. In addition, both the first and second order parameters of 
a quadratic fit to our entire sample are negative, again suggesting that the gradient 
steepens with increased galactocentric distance. However, we feel that given what 
appears to be natural scatter in the abundance data along with the presence of only a 
relatively small number of sample objects beyond about 17 kpc in the database, it is 
too early to conclude that the gradient is anything but linear with a single slope across 
the entire disk. 

Concerning the second point above, in order to achieve the fit described in point 1, we 
were forced to assume one-sigma abundance uncertainties that were 40% greater than the 
propagated ones; without this alteration, the value was too large. This result suggests 
the presence of real scatter in oxygen abundance among PNe at the same radial distance 
from the Galactic center. We suggested a few potential theories to explain this finding. 

Finally, we briefly considered the importance of achieving a consensus on the gradient 
value of the oxygen abundance distribution across the Galactic disk, as the current range 
in published slopes insufficiently constrains chemical evolution models and thus our under- 
standing of how the gradient originally formed and subsequently evolved. Related to this 
point is the need to clarify the actual form of the fitted model, e.g., linear, quadratic, discon- 
tinuous with multiple slopes, etc. The problem in doing this right now appears to be twofold. 
First, there is a paucity of oxygen measurements beyond 15 kpc from the Galactic center. 
Clearly, extending our sampling out farther will help to develop a consensus regarding the 
form of the abundance behavior, both qualitatively and quantitatively. But second, there is 
the problem of inherent, real scatter in the oxygen abundances whose cause is unknown but 
whose presence vastly complicates the study of the disk's oxygen distribution. This situation 
makes the study of the abundance gradient more difficult, of course, but it also brings to 
light a problem whose solution is bound to tell us something very interesting about how the 
level of interstellar oxygen in the Galaxy evolves. 

We wish to thank the staffs of the Apache Point Observatory and of KPNO and CTIO 
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also appreciate the help of Joe Rodgers and Eric Abraham regarding the methods for the 
statistical analysis of our data. RBCH, KBK, BB, and AJ also thank NSF for their gen- 
erous support through grants AST-0806577 to the University of Oklahoma, AST-080490 
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A. Abundance Comparison 



Our published work (including both the new object s reported here and o thers we have 
reported on elsewhere) includes 10 PNe in common with lWesson et al. I ( l2005l ). Our oxygen 
abunda nces for these 10 PNe agree with theirs on average to within 2%, at the worst differing 
by 17%. ICosta et al. I ( l2004l ) report on 11 objects in common with our total sample. Though 
their values are almost always higher than ours, agreement is reasonable for most of the 
objects, but for three (A12, Ml-13, and YC2-5) the difference exceeds a factor of two. 
Perhaps this is due to their having to rely on [O II] A7323 for the O"*" abundance, and not 
being able to measure A3727: for nine of these 11 objects (excluding for Ml-6, where we 
agree within 25%, and YC2-5, a high-ionization PNe with little 0+) the observed flux of 
A3727 is much stronger than that of A7323. 

Fig. |5] compares the oxygen abundances (0/H) for our entire published sample with 
those in SIO for the 126 objects in common. The gray line is the unity line. In general, 
the values lie within ±50% of the unity line; outliers are identified by name. In general, the 
points at low 0/H scatter around the unity line; at higher 0/H our values tend to be larger 
than SHlO's. Since SHIO do not explicitly list the abundance reference for each nebula, we 
cannot, in fact, trace the sources of the observed discrepancies. However, we can speculate on 
one possible cause, especially for highly-ionized nebulae: it could be the ionization-correction 
factor (IGF) for oxygen, which accounts for the presence of ions above O^^. Our IGF uses 
the (He"*" + He^"'")/He^ ratio with an exponent equal to 1, while a popular alternative, used 
by many groups, uses an exponent of 2/3. The latter IGF would reduce the final oxygen 
abundance compared with the former. In Fig. |5] the most egregious disparities occur for 
IG 418 and Gn 2-1, for which the SHIO abundance are much larger than ours. Both of these 
are low-excitation nebulae, so any IGF issue is moot. We do note that for two papers cited 
as abundance sources in SHIO, t he 0/H values in those papers are much closer to ours: 
Pottasch fc Bernard-Salas I ( 120061 ) derive O/H=3.5xl0~^ f o r IG 4 18, bringing it down much 



closer to the unity line; and for Gn 2-1, iGhiappini et al. I (|2009[ ) calculate 0/H= 4.7x10 
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which relocates it below and nearer to the unity line. 

As Fig. E] amply demonstrates, nebular abundances determined by different groups typ- 
ically exhibit discrepancies beyond their nominal uncertainties. These discrepancies presum- 
ably arise from choices made during observation, calibration and/or data analysis, each of 
which contains multiple alternatives that could all contribute to divergent final outcomes. 
Several of us (RBCH, KBK and BB) are embarking on a collaboration with a number of other 
nebular abundance groups to thoroughly investigate the sources of disagreements in derived 
abundances. As discussed in §4, these disparities currently represent a significant bottleneck 
hindering progress in using abundances to test ideas about galactic chemical evolution. 
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Table 1. Observing Log 



Nebula Name PNG Date Observed Total Exp (sec) Comments 



A 2 


122.1 -04.9 


Nov. 


A 8 


167.0 -00.9 


Nov. 


A 12 


204.0 -08.5 


Feb. 


A 14 


197.8 -03.3 


Feb. 


BV5-2 


121.6 -00.0 


Oct. 


BV5-3 


131.4 -05.4 


Oct. 


H3-75 


193.6 -09.5 


Feb. 


IC 289 


138.8 -1-02.8 


Feb. 


IC 1747 


130.2 -1-01.3 


Jan. 


IC 2149 


166.1 -1-10.4 


Jan. 


Kl-10 


229.6 -02.7 


Apr. 


K3-64 


151.4 -1-00.5 


Feb. 


K3-66 


167.4 -09.1 


Jan. 


K3-67 


165.5 -06.5 


Jan. 


K3-67 


165.5 -06.5 


Nov. 


K3-70 


184.6 -1-00.6 


Jan. 


K3-90 


126.3 -1-02.9 


Oct. 


K3-91 


129.5 -1-04.5 


Nov. 


K3-92 


130.4 -1-03.1 


Oct. 


K3-93 


132.4 -1-04.7 


Nov. 


K3-94 


142.1 -1-03.4 


Feb. 


K4-47 


149.0 -1-04.4 


Nov. 


K4-48 


201.7 -1-02.5 


Jan. 


Ml-4 


147.4 -02.3 


Jan. 


Ml-6 


211.2 -03.5 


Jan. 


Ml-7 


189.8 -1-07.7 


Jan. 


Ml-9 


212.0 -1-04.3 


Jan. 


Ml-11 


232.8 -04.7 


Feb. 


Ml-12 


235.3 -03.9 


Feb. 


Ml-14 


234.9 -01.4 


Apr. 


Ml-16 


226.7 -1-05.6 


Apr. 


M2-2 


147.8 -1-04.1 


Jan. 


M4-18 


146.7 -1-07.6 


Jan. 


NGC 1501 


144.5 -1-06.5 


Jan. 


NGC 2346 


215.6 -1-03.6 


Feb. 


PBl 


226.4 -03.7 


Apr. 


St3-1 


217.4 -1-02.0 


Feb. 


YC2-5 


240.3 -1-07.0 


Apr. 



2007 2700 

2007 3600 

2007 180 

2007 1800 PA=100° 

2007 3600 PA=45° 

2007 2700 

2007 1800 

2007 180 

2007 480 

2007 90 

2007 1800 

2007 1860 

2007 480 

2007 360 

2007 850 

2007 1800 

2007 2700 

2007 1500 

2007 2700 PA=45° 

2007 1800 

2007 960 

2007 1500 PA=40° 

2007 720 

2007 540 

2007 600 

2007 600 

2007 300 

2007 190 

2007 190 

2007 600 

2007 1200 

2007 900 

2007 420 

2007 240 

2007 180 

2007 900 

2007 240 

2007 600 



He2-15* 261.6 + 03.0 Nov. 2003 2100 B, 2400 R CTIO 

K3-61* 096.3 -I- 02.3 Aug. 2004 1200 B, 3000 R KPNO 

Ml-51* 020.9 - 01.1 Aug. 2004 1800 B, 2100 R KPNO 

NGC 6620* 005.9 - 06.2 Aug. 2004 2160 B, 4180 R CTIO 



*The line strength and abundance data for these non-APO objects have been included 
as an addendum to the online tables. 



Table 2. Fluxes and Intensities 



Line 


f(A) 




A 




A ft 


A 12 




A 1/1 


F(A)<= 


1(A) 


F(A) 


1(A) 


FfA) 




F(A) 


1(A) 


rn TTl A'?797 


0.292 


79.2 


116it26 


230 


416i93 


273 


382±86 


180 




[Ne III] A3869 


0.252 


76.6 


106±23 


60.8 


102±21 


82.3 


110±23 


38.1:: 


63.5±33.86:: 


He I + H8 A3889 


0.247 


13.8 


19.0±4.00 


12.5 


20.7±4.31 


13.9 


18.4±3.86 






n\.T^ TTTl \Qn£JO 

[INe ilij AoyDo 


0.225 


24.0^ 


DO n_Ln Qna 


23.8^ 


o'7 c;_Lin c;na 

0/ .o±lU.oy 


25.4=^ 


32.9±9.90=" 






He Ad97U 


0.224 


12.2^ 


16.3^ 


9.58^ 


15.1^ 


12. 5=^ 


16.2=" 






-tie 1 + xle 11 A4UZD 


0.209 


2.04:: 


z.Doihi.4z:: 














Me 11 A41UU 


0.188 


U.45z 


0.577^ 


r\ t A A 3. 

U.144 


0.211^ 


0.249=" 


0.309=" 






no A41U1 


0.188 


20.3^ 


zo.9±o.U4 


25.1^ 


36.8±7.U4 


19.5=" 


24.2±4.65=' 






Me 11 A4ooy 


U. Iz4 


U.oby 


1 noa 




n Qcna 


0.474=" 


0.546=" 


U.410 




TT- \/IO/irt 

M7 A4o4U 


0.124 


40.7^ 


47.8±8.3d 


32.6^ 


42.0±7.25 


42.0=" 


48.5±8.41=" 


47.2=" 


bU.7±lU.5U 


[S III] A9069 


-0.670 


92.2 


38.8±8.90 


102 


26.2±5.96 


38.9 


18.0±4.13 


22.6: 


5.84±2.12: 


P9 A9228 


-0.610 


9.70 


4.41±0.93 






4.78 


2.37±0.50 






[S III] A9532 


-0.632 


188 


83.0±18.01 


256 


71.1±15.30 


108 


52.4±11.38 


56.5: 


15.8±5.61: 


P8 A9546 


-0.633 


8.35 


3.68±0.80 






6.47 


3.12±0.68 






[C I] A9824 


-0.653 










11.9: 


5.62±2.03: 






[C I] A9850 


-0.654 










38.6 


18.2±4.08 






c 




0.56 




0.88 




0.50 




0.88 




Ha/H/3 




2.83 




2.86 




2.84 




2.86 




log Fh;3'= 




-13.59 




-14.21 




-12.42 




-14.82 





to 



="Deblended. 

'"ergs cm~^ in our extracted spectra 

"^The F(A) values have been scaled to H/3=100 using the observed Fh/3 values. Intensities of strong hnes have measurement uncertainties of 10%, 
single colons indicate uncertainties of 25%, double colons indicate uncertainties 50%, and triple colons indicate highly suspect measurements. 



Table 3. Ionic Abundances 



Ion 




A 2 




A R 

r\ 




A 12 




A 14 


used 


b 

Abundance 


used 


Abundance 


used 


Abundance 


T 

used 


Abundance 


He+ 


[O III] 


8.59±1.08(-2) 


[O III] 


7.23±3.66(-2) 


[O III] 


0.109±0.014 


[O III] 


0.172±0.087 


He+2 


[O III] 


4.00±1.31(-2) 


[O III] 


1.88±0.28(-2) 


[O III] 


2.14±0.32(-2) 


[0 III] 


2.40±1.25(-2) 


icf(He) 




1.00 




1.00 




1.00 




1.00 


00(6300) 


[N II] 




[N II] 




[N II] 


*2.58±0.67(-5) 


[N II] 


*2.96±1.24(-5) 




FN Til 

[IN lij 




[N TTl 

[IN lij 




[N TTl 

[IN ilj 




\N TTl 

[IN lij 




o° 


wm 




wm 




wm 


2.60±0.66(-5) 


wm 


3.31±1.37(-5) 


0+(3727) 


[N II] 


*3.43±0.75(-5) 


[N II] 


*9.92±7.11(-5) 


[N II] 


*7.40±1.89(-5) 


[N II] 


*5.35±2.44(-5) 


\j yi ozo ) 


["M TTl 


* 9 Q9_|_i fiof c:\ 


TNT TTl 
[IN lij 


*9 70-1- QO/' 


TNT TTl 


*fi 91 _|_9 7!^/ c;\ 


ri\T TTl 
[IN iiJ 




0+ 
w 


wm 


Q /II _|_n T'O/' 


wm 


Q 7n-l-R QS/' f^^l 

y. / uiCD.yot^-o ) 


wm 


7 '^9-1-1 S'?/' 


wm 


^ 9/1-1-9 A/l/ 


yj i^ouu / ) 


ro TTTl 


*9 K9-i-n /I ^ 


ro TTTl 


*9 oc;_i_n /in/' /1^ 
z.yOitU.4U^-4 J 


ro TTTl 


*9 99_|_n c:^;/ /i \ 
Z.ZZltU.OOl^-4 ) 


ro TTTl 


* 1 79_i_n 9Q/' /I ^ 

i. / ZltU.Zot^-4 J 




ro TTTl 


Z.OUltU.Ol(^-4j 


lr\ TTTl 

[U iiiJ 


Z. / yitU.oo(^-4 J 


\r\ TTTl 


i.UZitU.4U(^-4 j 


lr\ TTTl 


*l ^;c;_Ln T5/ /I ^ 
i.uoitU.Zo(-4J 




ro TTTl 


*9 fi9_|_n fifi/' /I \ 


ro TTTl 




rO TTTl 


*9 99_|_n c;f;/ a\ 
Z.ZZltU.ODl^-4 J 


fO TTTl 




n+2 


wm 


z.oyitU.Di(^-4j 


wm 


z.yiitU.oo(^-4 J 


wm 


1 Q-Un c; 1 /I \ 
Z.ioitU.Ol(^-4j 


wm 


1 7n-Ln 01 / A\ 
i. /UitU.Zi(-4J 




















s+ 


[N II] 


*5.31±0.84(-7) 


[N II] 


*9.23±3.54(-7) 


[N II] 


*8.85±1.89(-7) 


[N II] 


n.77±0.46(-6) 


S+(6716) 


[N II] 


5.31±0.85(-7) 


[N II] 


9.25±3.53(-7) 


[N II] 


8.85±1.89(-7) 


[N II] 


1.77±0.46(-6) 


S+(6731) 


[N II] 


5.31±0.84(-7) 


[N II] 


9.21±3.56(-7) 


[N II] 


8.85±1.89(-7) 


[N II] 


1.76±0.47(-6) 


S+2(9069) 


[S III] 


*4.23±1.55(-6) 


[O III] 


*4.01±0.86(-6) 


[S III] 


*2.10±0.65(-6) 


[O III] 


*8.92±3.17(-7) 


3+2(6312) 


[S III] 


*4.23±1.55(-6) 


[O III] 




[S III] 


*2.10±0.65(-6) 


[O III] 




g+2 


wm 


4.23±1.55(-6) 


wm 




wm 


2.10±0.65(-6) 


wm 




icf(S) 




1.37±0.06 




1.19±0.08 




1.19±0.03 




1.19±0.05 



^ wm indicates the calculated mean value for the ionic abundance weighted by the observed flux 
* indicates those values included in the weighted mean 



Table 4. Temperatures and Densities 



Parameter 


A 2 




A 8 




A 12 


A 14 




Value"" 


Notes'^ 


Value 


Notes 


Value Notes 


Value 


Notes 


T[o III] 


11710±649 




10000 


Default. 


11430±618 


10000 


Default. 


T[N II] 


10300 


Default. 


11580±2791 




12160±932 


13270±2306 




T[o II] 


9453±3083 




6425±1409 




10850±2680 


11450±2833 




T[s III] 


11050±1923 


Used 9069. 






12460±1275 Used 9532. 






Ne,[S II] 


73.2±179.30 




75.2±169.90 




301±244 


55.1±164.50 





^Temperatures and densities given in kelvin and cm ^ respectively. 

^If [N II] A5755 is unavailable, Tj^ jj] is estimated from Tjq m] (Kwitter &: Henry 2001). Since telluric absorption tends 
to affect only one or the other of [S III] A9069 and A9532, if the observed flux ratio A9532/A9069 > 2.48 (the theoretical 
value) then A9532 is used in the Tjs m] calculation; otherwise A9069 is used. Default T[n h] and ^c[c\ iii] values and the 
high density limit for Ne[s ii] are based on criteria discussed in Kwitter &; Henry (2001). 
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Table 5. Total Elemental Abundances 



Parameter ^ 


A 2 


A 8 


A 12 


A 14 


Solar Ref 


Orion Ref ^ 


He/H 


0.126±0.018 


9.11±3.68(-2) 


0.130±0.015 


0.196±0.088 


8.51(-2) 


9.73(-2) 


N/H 


1.34±0.47(-4) 


1.48±0.40(-4) 


1.32±0.37(-4) 


1.05±0.24(-3) 


6.76(-5) 


5.37(-5) 


N/0 


0.313±0.088 


0.302±0.126 


0.380±0.090 


4.12±1.27 


0.138 


0.100 


0/H 


4.29±1.05(-4) 


4.89±1.19(-4) 


3.48±0.71(-4) 


2.54±0.47(-4) 


4.90(-4) 


5.37(-4) 


Ne/H 


9.93±2.71(-5) 


1.76±0.56(-4) 


1.08±0.25(-4) 


9.79±5.48(-5) 


8.51(-5) 


1.12(-4)) 


Ne/0 


0.231±0.037 


0.361±0.063 


0.311±0.051 


0.385±0.201 


0.174 


0.209 


S/H 


6.52±2.37(-6) 


5.88±1.18(-6) 


3.54±0.97(-6) 


3.15±0.69(-6) 


1.32(-5) 


1.66(-5) 


S/0 


1.52±0.61(-2) 


1.20±0.41(-2) 


1.02±0.31(-2) 


1.24±0.33(-2) 


2.69(-2) 


3.09(-2) 


Cl/H 






4.14±1.58(-8) 




3.16(-7) 


2.88(-7) 


Cl/0 






1.19±0.48(-4) 




6.46(-4) 


5.37(-4) 


Ar/H 


1.98±0.43(-6) 


2.66±0.71(-6) 


1.30±0.26(-6) 


1.46±0.28(-6) 


2.51(-6) 


4.17(-6) 


Ar/0 


4.61±0.96(-3) 


5.44±1.47(-3) 


3.72±0.78(-3) 


5.75±1.22(-3) 


5.13(-3) 


7.76(-3) 



'^all abundances are in linear form 

^Solar values taken from Asplund et al. 2009. 

'^Orion values taken from lEsteban et al. I (j2004l ) 
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Table 6. Data Used In Analysis 



PN Name Dgsv' Rssv^ 12+log(0/H) 



A 1 9 


9 04 


1 n 44 


9 Dfi 


1 n 4(S 


R ^4 


rn9 1 


o.ou 


4 (S7 


8 nn 

o.uu 


n 7(S 


R sn 

o.ou 


r,ni, 1 


o.oo 


(\ 14 


(S n4 




8 fin 

o.uu 


Fp-1 
r g± 


3 1 


7 96 


3 1 3 


7 96 


O. O X 


H9 1 8 




4 (S^ 






R 7^ 
O. 1 o 




9 70 


1114 

± ± . ±4: 






R ^R 
o.oo 




9 OR 


4"^ 


^ 04 


4R 


R R9 

O .OZ( 


IIUU 


1 


u. ou 






o.uu 


Hp9-1 1 1 

LLKJZj J- J- J- 


O .t:0 


U. UZ( 


^1 


8 ^1 


R AA 

O .tit: 


Hp9 11"^ 


1 

± .c/O 


7 HQ 






R ^n 

o.ou 


Hp9 141 


O .tit: 


nn 


47 


^ QR 


R RO 
o.ou 


Hp9 1 


9 14 


Q 

y. uu 


9 1 7 

Z,. 1 / 


Q n7 


R "^7 

O.O 1 


Hp9 1 '^R 


1 Q 74 


1 9"^ 


1 Q QD 


1 "^R 
±o .oo 


R 4"^ 

O .4:0 


Hp9 91 


7 


1 n (^7 


1 9 


14 49 


R "^R 
o.oo 


Hp9 '^T 


9 QR 

Zi . c/O 


R 7R 

O. 1 o 




R 7R 

O. 1 o 


O .C/O 


Hp9 dR 


1 n 


o. oo 






O.OO 


Frp9 f^f^ 


9 f^n 

z.uu 


o. lu 






R 7f> 
O. 1 u 


Hul-z 


1.48 


8.54 






8.23 


Hu2-1 


2.5 


7.22 






8.43 


IC1297 


3.76 


5.01 


4.91 


3.94 


8.80 


IC1747 


2.94 


10.64 


2.99 


10.68 


8.57 


IC2149 


1.59 


10.03 


3.26 


11.64 


8.35 


IC2165 


2.21 


10.22 


3.65 


11.42 


8.40 


IC2501 


2.62 


8.40 






8.57 


IC2621 


1.87 


8.00 


4.87 


8.09 


8.61 


IC3568 


2.71 


9.91 


2.74 


9.93 


8.51 


IC418 


0.61 


8.96 






8.28 


IC4593 


3.19 


6.40 


3.23 


6.38 


8.62 


IC4776 


3.9 


4.71 


4.97 


3.67 


8.61 


IC5217 


4.65 


10.40 


5.37 


10.84 


8.51 
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Table 6 — Continued 



PNName Dgsv' Rssv^ 12+log(0/H) 



J320 


6.06 


14.22 


6.12 


14.27 


8.32 


J900 


2.76 


11.19 


4.90 


13.30 


8.47 


K3-61 


7.33 


11.81 


7.41 


11.87 


8.60 


K3-64 


9.82 


17.75 






8.42 


K3-67 


3.4 


11.80 


3.43 


11.83 


8.12 


K3-70 


12.12 


20.60 


19.40 


27.88 


8.08 


K3-91 


8.44 


15.30 


8.52 


15.37 


8.34 


K3-92 


6.81 


13.91 


6.87 


13.96 


8.50 


K3-93 


7.45 


14.58 






8.40 


K3-94 


6.51 


14.20 


6.57 


14.26 


8.51 


K4-47 


8.53 


16.39 


8.61 


16.46 


7.69 


K4-48 


7.77 


15.97 


13.50 


21.61 


8.51 


Ml-13 


5.32 


12.48 


5.38 


12.53 


8.70 


Ml-14 


3.95 


11.24 


6.42 


13.27 


8.46 


Ml-16 


5.45 


12.84 


9.77 


16.74 


8.52 


Ml-17 


7.36 


14.42 


10.80 


17.56 


8.69 


Ml-34 


4.99 


3.54 


5.05 


3.48 


8.70 


Ml-4 


2.99 


11.13 


6.60 


14.50 


8.40 


Ml-40 


1.81 


6.71 






8.75 


Ml-42 


5.46 


3.08 


5.51 


3.03 


8.55 


Ml-5 


2.92 


11.41 






8.16 


Ml-50 


4.36 


4.43 


5.98 


3.12 


8.73 


Ml-51 


1.75 


6.89 


2.29 


6.41 


8.94 


Ml-54 


3.78 


4.99 


3.81 


4.96 


8.64 


Ml-7 


5.91 


14.31 


5.96 


14.36 


8.67 


Ml-74 


4.12 


6.81 






8.62 


Ml-8 


3.39 


11.55 


3.42 


11.58 


8.65 


Ml-80 


5.5 


11.44 


5.55 


11.47 


8.87 


Ml-9 


4.88 


12.89 


11.60 


19.31 


8.28 


M2-52 


4.41 


10.46 


4.45 


10.49 


8.33 


M2-55 


2.21 


9.67 


2.23 


9.68 


8.93 
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Table 6 — Continued 



PN Name Dgsv" Rssv^ 12+log(0/H) 



M3-15 


0.29 


4.47 








8.87 


M3-2 




8.82 


14.91 


8, 


.93 


15.01 


7.83 


M3-3 




5.75 


13.33 


5, 


.81 


13.38 


8.47 


M3-4 




6.33 


12.82 


6.39 


12.87 


8.58 


M3-5 




7.14 


13.18 


7.21 


13.24 


8.41 


M3-6 




3.29 


9.92 


4, 


.22 


10.47 


8.70 


Mel-1 


4.62 


6.77 


6, 


.76 


6.93 


8.69 


Me2-2 


5 


10.55 








8.41 


Mz2 




2.34 


6.60 


2, 


.36 


6.58 


8.95 


Mz3 




1.27 


7.41 


1, 


.45 


7.26 


8.16 


NGC 


2346 


1.36 


9.64 


1, 


.37 


9.64 


8.60 


NGC 


2371 


1.54 


9.93 


1, 


.55 


9.94 


8.79 


NGC 


2392 


1.25 


9.64 


1, 


.26 


9.65 


8.39 


NGC 


2438 


1.2 


9.29 


1, 


.22 


9.30 


8.61 


NGC 


2440 


1.35 


9.34 


1, 


.36 


9.35 


8.62 


NGC 


2452 


2.81 


10.08 


2, 


.84 


10.10 


8.87 


NGC 


2792 


3.02 


9.23 


3 


.05 


9.24 


8.82 


NGC 


2867 


1.84 


8.44 


2, 


.23 


8.48 


8.64 


NGC 


3195 


1.96 


7.85 


1, 


.98 


7.85 


8.71 


NGC 


3211 


2.87 


8.17 


2, 


.90 


8.17 


8.84 


NGC 


3242 


1.08 


8.69 


1, 


.09 


8.69 


8.57 


NGC 


3587 


0.62 


8.79 


0, 


.62 


8.79 


8.54 


NGC 


3918 


1.01 


8.13 


1, 


.64 


7.96 


8.67 


NGC 


5307 


3.2 


6.79 


3, 


.24 


6.78 


8.50 


NGC 


5315 


1.24 


7.78 








8.64 


NGC 


5882 


1.68 


7.15 


2, 


.36 


6.65 


8.68 


NGC 


6210 


2.03 


7.41 


2, 


.28 


7.29 


8.68 


NGC 


6302 


0.53 


7.98 


0, 


.74 


7.77 


8.21 


NGC 


6309 


2.53 


6.10 


2, 


.74 


5.90 


8.77 


NGC 


6369 


0.66 


7.84 


1, 


.09 


7.42 


8.71 


NGC 


6445 


1.37 


7.15 


1, 


.38 


7.14 


8.83 
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Table 6 — Continued 



PN Name Dgsv" Rssv^ 12+log(0/H) 



NGC 


6439 


4.12 


4.54 


6.33 


2.61 


8.89 


NGC 


650 


0.74 


8.99 


7.46 


14.41 


8.80 


NGC 


6537 


0.9 


7.62 






8.32 


NGC 


6563 


1.63 


6.88 


1.65 


6.86 


8.66 


NGC 


6565 


4.62 


3.91 


4.66 


3.87 


8.75 


NGC 


6567 


2.37 


6.20 


3.61 


5.02 


8.36 


NGC 


6572 


0.71 


7.94 


1.74 


7.16 


8.58 


NGC 


6578 


2.31 


6.25 


3.64 


4.97 


8.83 


NGC 


6620 


8.75 


0.89 


8.84 


0.92 


8.86 


NGC 


6629 


1.95 


8.50 


2.37 


6.18 


8.78 


NGC 


6720 


0.87 


6.89 


0.88 


8.15 


8.79 


NGC 


6741 


2.05 


6.40 


3.73 


5.79 


8.61 


NGC 


6751 


2.56 


6.20 


2.89 


6.15 


8.63 


NGC 


6790 


1.54 


7.35 






8.60 


NGC 


6803 


2.99 


6.80 


5.27 


6.19 


8.74 


NGC 


6826 


1.58 


8.47 


1.59 


8.47 


8.57 


NGC 


6853 


0.26 


8.50 


0.26 


8.38 


8.72 


NGC 


6881 


2.47 


8.19 


5.34 


8.75 


8.65 


NGC 


6884 


2.11 


8.47 


3.83 


8.82 


8.69 


NGC 


6886 


3.1 


7.46 


4.35 


7.37 


8.65 


NGC 


6891 


3.19 


7.13 


3.61 


7.04 


8.57 


NGC 


7008 


0.86 


8.59 


0.87 


8.59 


8.76 


NGC 


7009 


1.2 


7.74 


1.33 


7.66 


8.72 


NGC 


7026 


1.9 


8.68 


2.35 


8.78 


8.81 


NGC 


7027 


0.27 


8.48 






8.52 


NGC 


7293 


0.16 


8.40 


0.16 


8.43 


8.66 


NGC 


7354 


1.27 


8.97 


1.70 


9.16 


8.71 


PB6 




4.38 


8.94 


4.42 


8.95 


8.69 


PC14 




5.74 


4.01 


5.80 


3.98 


8.84 


Pel-1 


8 


3.56 


5.58 






8.57 


St3-1 




5.80 


13.57 






8.58 
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Table 6 — Continued 



PN Name 






Dssv'^ 


R-ssv*^ 12- 


hlog(0/H) 


Th2-A 


2.45 


7.32 


2.47 


7.31 


8.71 



'^Most of the heliocentric distances in this column 



were taken directly from ICahn et al. I (119921 ). The 



sources for the remaining objects are provided in Sec- 
tion 4. 

*^Galactocentric distances were computed from the 
heliocentric distances in column 2 using the relation: 
R = {R| + [cos(b) X D]2 - 2 X Rq X D X cos(l) x 
cos(b)}^/^, where R and D are galactocentric and 
heliocentric distances, respectively, and b and 1 are 
galactic latitude and longitude, respectively. We as- 
sumed that Rq = 8.5 kpc. 



'^Heliocentr ic distances in this column were taken 
directly from IStanghellini et al. I ( l2008l ). 
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Table 7. Trial Least Squares Fits 



f(0) 


f(Rg) 


b 




1 





-0.041 


1.95E-28 


1 


0.1 


-0.051 


4.70E-16 


1.25 





-0.041 


3.77E-10 


1.25 


0.1 


-0.049 


4.92E-06 


1 


0.2 


-0.066 


0.00071 


i.o 





-U.04i 


O.OOoO 


1.5 


0.1 


-0.047 


0.050 


1.25 


0.2 


-0.061 


0.13 


1.75 





-0.041 


0.42 


1.75 


0.1 


-0.045 


0.71 


1.5 


0.2 


-0.056 


0.75 


1 


0.3 


-0.074 


0.78 


2 





-0.041 


0.96 


1.25 


0.3 


-0.070 


0.98 


2 


0.1 


-0.045 


0.99 


1.75 


0.2 


-0.053 


0.99 


2 


0.2 


-0.051 


1 


1.5 


0.3 


-0.066 


1 


2 


0.3 


-0.059 


1 


1.75 


0.3 


-0.062 


1 



Table 8. Oxygen Gradients 



Sample'^ 


n 


a 


b 








r 


Pr 


Total 


124 


9.09±.05 


-0.058±.006 


121.8 


1.00 


0.49 


-0.54 


3.3(-ll) 


Type I 


48 


9.10±.07 


-0.061±.008 


62.7 


1.36 


0.051 


-0.59 


3.5(-6) 


Type II 


75 


9.05±.08 


-0.053±.010 


58.6 


0.80 


0.89 


-0.44 


4.8(-5) 


R< lOkpc 


85 


9.04±.09 


-0.054±.013 


73.3 


0.88 


0.77 


-0.21 


4.8(-2) 


R> lOkpc 


38 


9.92±10.10 


-0.12±.14 


24.2 


0.67 


0.93 


-0.39 


1.2(-2) 


H II Regions 


73 


9.24±.07 


-0.078±.007 


63.5 


0.89 


0.72 


-0.80 


1.2(-20) 


B V Stars 


24 


9.63±.ll 


-0.080±.01 


33.3 


1.51 


0.058 


-0.64 


2.6(-4) 



^The first five groups refer to the PN sample presented here, where the abundance 
uncertainty for each object was increased by 40% and the error in the galactocentric 
di stance was assumed to b e 20%. The H II region sam ple w as assembled from studie s 



bv lAffierbach et al. I (|l997l ). lvnchez fc Esteban I (119961 ^. an dlPeharveng et al. I(l2000l ). 



and the B V st a r sarn ple was compiled from papers by ISmartt et al. 
Rolleston et"al~l (|200ol ). 



(2001 



) and 



Table 9. Oxygen Gradient Comparison 



Source^ N Alog(0/H)/AR (dex/kpc) Object Type ARg (kpc) 



This Paper 


124 


-0.058±.006 


PN 


0.9-21 


SIO 


145 


-0.023±.006 


PN 


3-21 


H04 


79 


-0.037±.008 


PN 


2-17 


C04 


80 


-0.05 ± • • • 


PN 


4-14 


P06 


83 


-0.016±.008 


PN 


0.5-15 


R06 


70 


-0.060±.010 


H II (optical) 


5-18 


R06 


68 


-0.041±.014 


H II (FIR) 


0.1-15 


P09 


265 


-0.051±.004'^ 


Cepheids 


5-17 


SR97 


47 


-0.07±.01 


B V Stars 


6-18 



C04: 
R06: 



SIO: IStanghellini Haywood 



Costa et al. 



Rudolph et al 



Smartt &: Rolleston 



mm- 



120061 1: 
(|l997f ) 



(2010 


); H04: 


Henrv et al. 


9: 


Perinotto k, Morbidelli 


Pedicelli et al. ( 


2009 





(2004); 
(|2006l l: 



SR97: 



'^Alog(Fe/H)/AR (dex/kpc) 
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R (kpc) 



Fig. 1. — 12+log(0/H) versus galactocentric distance, R, in kpc for PNe in our total sample. 
Least squares fits are indicated by line type for the total sample (solid), objects for R<10kpc 
(dashed), R>10kpc (dot-dashed), while a quadratic fit is shown with a dotted line. 
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Fig. 2. — 12+log(0/H) versus galactocentric distance, R, in kpc for objects in our sample. 
Type I PNe are symbolically distinguished from Type lis as shown in the legend. Least 
squares fits for Type I and Type II are shown with dot-dashed and solid lines, respectively. 
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Fig. 3. — 12+log(0/H) versus galactocentric distance, R, in kpc. PNe are separated accord- 
ing to whether they are greater than or less than 300 kpc from the Galactic plane in the 
vertical direction. 
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• PNe 

O H II Regions 

X B V Star 

■ ■ B V Stars 

PNe 

H II Regions 




15 

R (kpc) 



30 



Fig. 4. — 12+log(0/H) versus galactocentric distance, R, in kpc. Planetary nebulae are 
indicated with filled circles, while H II regions are indicated with open circles, and stars are 
indicated with stars. The lines are least squares fits to the respective data types, as indicated 
in the legend. The sources for the samples of PNe, H II regions, and B V stars are the same 
as those provided in Table [81 
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2 4 6 8 10 



O/H (This Paper) x 10* 



Fig. 5. — O/H X 10^ abundances from lStanghellini fc Haywood I (120101 . vertical axis) and the 
present sample (horizontal axis) for objects common to both. The 1:1 track is shown with a 
straight line. 



